Non equilibrium dynamics of disordered systems : understanding the broad 
continuum of relevant time scales via a strong-disorder RG in configuration space. 
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We show that an appropriate description of the non-equilibrium dynamics of disordered systems 
is obtained through a strong disorder renormalization procedure in configuration space, that we 
define for any master equation with transitions rates W {C ^ C') between configurations. The idea 
is to eliminate iteratively the configuration with the highest exit rate Wout{C) — X^c' (C — > C) 
to obtain renormalized transition rates between the remaining configurations. The multiplicative 
structure of the new generated transition rates suggests that, for a very broad class of disordered 
systems, the distribution of renormalized exit barriers defined as BoutiC) = — In WoutiC) will become 
broader and broader upon iteration, so that the strong disorder renormalization procedure should 
become asymptotically exact at large time scales. We have checked numerically this scenario for the 
non-equilibrium dynamics of a directed polymer in a two dimensional random medium. 



The non-equiUbrium dynamics of disordered systems 
usually displays a broad continuum of relevant time 
scales, that give rise to a lot of striking properties such as 
aging, rejuvenation and memory that have been studied 
a lot both experimentally and theoretically (see [l| and 
references therein). In finite dimensions, these effects can 
be understood via the growth of some coherence length 
l{t) that separates the smaller lengths I < l{t) which are 
quasi-equilibrated from the bigger lengths I > l{t) which 
are completely out of equilibrium. The slow nature of 
the dynamics then reflects the fact that equilibration on 
larger length scales requires to overcome larger and larger 
barriers. Within the droplet scaling theory proposed 
both for spin-glasses 0, Q and for the directed polymer 
in a random medium the non-equilibrium dynamics 
is activated with barriers scaling as power law B{1) ~ 
of the length scale I. The typical time associated to scale 
I then grows as an exponential h\ttyp{l) ~ B{1) ~ , or 
equivalently, the characteristic length-scale l{t) associ- 
ated to time t grows only logarithmically l{t) ~ (Int)"'^^'''. 
In numerical studies, this logarithmic behavior has re- 
mained controversial because the maximal equilibrated 
length Imax measured at the end of the simulations is 
usually rather small, so that many fits of the data are 
possible. For instance, in Monte-Carlo simulations of 
2D or 3D random ferromagnets [S] or spin-glasses Q, 
the maximal equilibrated size is usually only of order 
Imax 10. A noteworthy exception is the simulation of 
an elastic line in a random medium where sizes of order 
hnax ^ 100 have been measured 0] with the conclusion 
that the length l{t) grows logarithmically with a barrier 
exponent tp ^ 0.49 ( whereas power- law fits l{t) ~ t^/^ 
are excluded at large times). Note that the difficulties 
met in dynamical Monte Carlo simulations of disordered 
systems comes precisely from the presence of a contin- 
uum of relevant time scales ranging from the microscopic 
scale of single moves to the equilibrium time of the full 
system. Then even faster-than-the-clock Monte Carlo al- 
gorithms @], where each iteration leads to a movement, 
become inefficient because they face the 'futility' prob- 
lem : the number of different configurations visited 



during the simulation remains very small with respect to 
the accepted moves. The reason is that the system vis- 
its over and over again the same configurations within a 
given valley before it is able to escape towards another 
valley where it will be trapped even longer! 

In this paper, we argue that the appropriate descrip- 
tion of these dynamics with a broad continuum of rele- 
vant time scales requires some renormalization where the 
smaller time scales are successively integrated out to ob- 
tain the properties of the large-time dynamics. Moreover, 
since we expect that in these systems, disorder dominates 
over thermal fluctuations at large scales, the most ap- 
propriate renormalization scheme is a so called 'strong 
disorder renormalization (RG) procedure' (see [13] for a 
review). This very specific type of RG, which was in- 
troduced by Ma and Dasgupta [ll| and developed by 
D.S. Fisher [l2| in the field of quantum spin chains, has 
been then successfully applied to various classical disor- 
dered dynamical models, such as random walks in ran- 
dom media p^ . reaction-diffusion in a random medium 
p^ . coarsen ing dynamics of classical spin chains 15 1, 
trap models [l^|, absorbing state phase transitions 17| . 
zero range processes [Hj, exclusion processes [l^. In all 
these cases, the strong disorder RG rules have been for- 
mulated in real space, with specific rules depending on 
the problem. In this paper, we show that for more com- 
plex systems where the formulation of strong disorder 
RG rules has not been possible in real space, it is never- 
theless possible to formulate strong disorder RG rules in 
configuration space. Moreover, this formulation in con- 
figuration space is very general since it can be defined 
for any master equation describing the evolution of the 
probability -Pt(C) to be in configuration C at time t 
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The notation W (C — > C) represents the transition rate 
per unit time from configuration C to C, and 



Wo. 



(2) 
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represents the total exit rate out of configuration C. 
The two important properties of this master equation 
are the following (i) the exit time r from configuration 
C is distributed with the exponential law Pq"*{t) — 
Wout (C) e"'^'^""*^''^ (ii) the new configuration C where 
the jumps jumps at time r when it leaves the configura- 
tion C is chosen with the probability ttc (C) = 

For dynamical models, the aim of any renormalization 
procedure is to integrate over 'fast ' processes to obtain 
effective properties of 'slow' processes. The general idea 
of 'strong renormalization' consists in eliminating iter- 
atively the 'fastest' process. For the master equation of 
Eq. [1] we thus define the strong disorder renormalization 
in configuration space by the iterative elimination of the 
configuration with the highest exit rate (Eq. [2]). Let us 
call this configuration C* , and its exit rate W*^f. 



Wont (C*) = niaxc [Wont (C)] 



(3) 



We now have to compute the 'new' effective transitions 
rates M^"'^"'(C C) among the remaining configura- 
tions in terms of the 'old' transitions rates W°''^{C — > C) 
where the decimated configuration C* was still present. 
The only changes occur for the configurations called 
here (Ci,C2, ...,C„) that were related via positive rates 
W°^'^{C* Ci) > and W^°''^(C, ^ C*) > to the dec- 
imated configuration C* (here we will assume, for the 
simplicity of the discussion, and because it is usually the 
case in statistical physics models, that if a transition has 
a strictly positive rate, the reverse transition has also a 
strictly positive rate; but of course the renormalization 
rules can be simply extended to other cases). The 2n 
rates W°^'^{C* C,) and W°^'^{C^ C*) with i = 1, ..,n 
have to be eliminated, after taking into account their 
effects on transitions between pairs of neighbors of C*. 
For each neighbor configuration Ci with i e (l,..,n), 
the renormalized rate to go to the configuration Cj with 
j € (1, .., n) and j =/= i reads 



W^"'="'(C, ^Cj) = 
+W''^'^{Ci C*) 



Wl 



(4) 



The first term represents the 'old' transition rate (possi- 
bly zero), whereas the second term represents the tran- 
sition via the decimated configuration C* : the factor 
W°^'^{Ci — > C*) takes into account the transition rate to 
C*, whereas the second factor represents the probability 
to make a transition towards Cj when in C* . Note that 
the rule of Eq. 2] has been recently proposed in (2^ to 
eliminate 'fast states' from various dynamical problems 
with two very separated time scales. The physical inter- 
pretation of this rule is as follows : the time spent in the 
decimated configuration C* is neglected with respects to 
the other time scales remaining in the system. The va- 
lidity of this approximation within the present renormal- 
ization procedure will be discussed in detail below. To 
finish the decimation of C*, we have now to update the 



exit rates out of the neighboring configurations Ci 
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Since the only changes come from the rate towards C* 
that has disappeared and from the rates towards j G 
(1, ..,n) with j ^ i that have changed according to Eq. 
m one obtains 



w:Li{c.) 



W'^iC -> c* 



^rold^Q* _> Ci) 
w* t 



The physical meaning of this rule is the following. The 
exit rate out of the configuration d decays because the 
previous transition towards C* could lead to an immedi- 

ate return towards Ci with probability ^ —. Af- 
ter the decimation of the configuration C*, this process is 
not considered as an 'exit' process anymore, but as a res- 
idence process in the configuration Ci. This point is very 
important to understand the meaning of the RG proce- 
dure : the remaining configurations at a given stage are 
'formally' microscopic configurations of the initial master 
equation, but each of these remaining microscopic config- 
uration actually represents some 'valley' in configuration 
space that takes into account all the previously decimated 
configurations. Note that in practice, the renormalized 
rates W{C — > C) can rapidly become very small as a 
consequence of the multiplicative structure of the renor- 
malization rule of Eq[51 So the appropriate variables are 
the logarithms of the transition rates, called 'barriers' 
from now on. The barrier B{C — + C) from C to C is 
defined by B(C C) = -lnM^(C ^ C) and similarly 
the exit barrier out of configuration C is defined by 



Bout{C) = -\nWout{C) 



(6) 



As mentioned above, the approximation made in the 
renormalization rule of Eq. 3] consists in neglecting the 
time spent in the decimated configuration C* with re- 
spects to the other time scales remaining in the system. 
In our present framework, this means that the maximal 
exit rate chosen in Eq [3] should be well separated from 
the exit rates of the neighboring configurations Cj. The 
crucial idea of 'infinite disorder fixed point' [Si is 
that even if this approximation is not perfect during the 
first steps of the renormalization, this approximation will 
become better and better at large RG scale if the prob- 
ability distribution of the remaining exit rates becomes 
broader and broader upon iteration. More precisely, if 
the renormalization scale F is defined as the exit barrier 
of the last eliminated configuration F = Bout{C*), one 
expects that the probability distribution of the remain- 
ing exit barriers Bout > F will converge towards some 
scaling form 



PriBout-T) 



1 p f Bout 



a(F) 



a(F) 



(7) 



where P is the fixed point probability distribution, and 
where cr(F) is the appropriate scaling factor. The notion 
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of 'infinite disorder fixed point' means that the width 
<t(T) grows to infinity with the RG scale F. For instance, 
in previously known cases of infinite disorder fixed points 
where calculations can be done explicitly [fO], the scale 
ct(F) grows linearly a(T) ^ F, and the fixed point distri- 
bution is an exponential P{x) = e^^. Whenever the flow 
is towards an 'infinite disorder fixed point', the strong 
disorder renormalization procedure becomes asymptoti- 
cally exact at large RG scales. For our present problem, 
the convergence towards an 'infinite disorder fixed point' 
will depend on the initial condition of the transition rates, 
i.e. on the model ( and on the temperature if there are 
phase transitions). However, the form of the renormal- 
ization rule of Eq H] is sufficiently similar to the usual 
Ma-Dasgupta rules [l^] to think that the convergence 
towards some infinite disorder fixed point should be real- 
ized in a very broad class of disordered systems in their 
glassy phase. In practice, it can be checked numerically 
for each model of interest. 
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consider the usual Metropolis dynamics at temperature 
T — 0.5 defined by the transition rates 

W{C^ C) = S^cfi'> min (l, e-(^^(c')-s(C))/T^ 

where the factor (5<c,c'> means that the two configu- 
rations are related via the move of a single monomer 
hx hx i2. So initially each configuration has at most 
L neighbors corresponding to single moves. However dur- 
ing the renormalization, many new transition rates will 
be generated via Eq. [H as in real-space strong disor- 
der RG studies of quantum models in dimension d > 1 
[21L [2^ . Here, to validate our approach, we have decided 
to follow exactly the full RG flow, without disregarding 
any new transition rate. As a consequence, we have been 
able to study numerically only moderate lengths L < 9 
( 2^ < 512 configurations ) with a sufficient statistics of 
Us > 10^ disordered samples (we have data up to i = 11 
with rig = 500 samples, but histograms are too noisy). 
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FIG. 1: (Color online) Flow of the probability distribution 
Pr{Bout — r) of the renormalized exit barriers (see Eq. [7]) 
as the RG scale grows F = 3, 4, 5, 6, 7. Inset : growth of the 
width o-(F) with the RG scale F. (Data obtained for a di- 
rected polymer of length L = 9 in a two-dimensional random 
medium with 2^ — 512 configurations, the statistics is over 
Us = 15.10'* disordered samples). 

As an example of application, we consider the di- 
rected polymer in a two-dimensional random medium, 
a model first introduced to describe interfaces in ran- 
dom ferromagnets [2^ (see [l^l for a review). The stat- 
ics is well described by the Fisher-Huse droplet the- 
ory [3| as checked by detailed numerical studies [1, . 
We use the following discrete formulation for a poly- 
mer of length L attached at the origin : the 2^ con- 
figurations (hi, ...hi,) correspond to the random walks 
hx — hx-i = ±1 starting at /iq = 0. The energy of a con- 
figuration IS E = X]k=i ^(■^' where the site random 
energies e{x,h) are Gaussian p{t) = e~'^^/^/V27r. We 



FIG. 2: (Color online) Statistics of the equilibrium time over 
the samples for a directed polymer in a two-dimensional ran- 
dom medium : probability distribution QL^Teq = In teg) for 
length L = 5,6,7,8,9 (see Eq. |8]). Inset : rescaled distribu- 
tion Q{u) in log representation to see the tails. 

We find (see Fig 1) that the probability distribution 
of renormalized exit barrier of Eq. [7] flows towards an 
"inflnite disorder" fixed point, with a width growing as 
(t(F) ~ F and a rescaled probability which is extremely 
close to the exponential P{x) ~ e~^. Note that this 
type of renormalized distribution seems extremely robust 
within strong disorder RG since they hold for exactly in 
soluble models in c? = 1 [l3| and have been also found 
numerically in quantum models in dimension d > 1 [2l|. 
In a finite sample, the typical equilibrium time t^q can 
be obtained from the last decimated barrier T^q leading 
to a single surviving configuration via Fg^ = Intgq. Its 
probability distribution Qii^eq = Inigg) over the disor- 
dered samples of size L is shown on Fig. 2 for various L. 
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The convergence towards a fixed rescaled distribution 



leys above a given life-time, as well as aging properties. 



QLi^eq 



A{L) 



Q\u = 



ML) 



(8) 



is rapid (sec inset of Fig. 2) but the sizes studied are not 
sufhcient to obtain, via the average TeqiL) ~ or the 
width A{L) ~ L"^ , a reliable measure of the asymptotic 
barrier exponent ^, whose value has remained contro- 
versial (see [26| for a recent summary). We hope in the 
future to propose simplified ways of following the impor- 
tant rates of the renormalized flow to reach bigger system 
sizes f??! . We also intend to study other properties of the 
RG flow, in particular the structure of the evolving set of 
remaining configurations that label the metastable val- 



In conclusion, we have proposed to describe the 
non-equilibrium dynamics of disordered systems via a 
strong disorder renormalization procedure in configura- 
tion space, that we have defined for any master equa- 
tion. We have argued that for a very broad class of 
disordered systems, the distribution of renormalized exit 
barriers should become broader and broader upon itera- 
tion, so that the strong disorder renormalization proce- 
dure should become asymptotically exact at large time 
scales. We have checked this scenario numerically for the 
directed polymer in a two-dimensional random medium 
and we intend to study other disordered models in the 
future [27|. 
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